Double layer in ionic liquids: Overscreening vs. crowding 
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We develop a simple Landau- Ginzburg- type continuum theory of solvent-free ionic liquids and 
use it to predict the structure of the electrical double layer. The model captures overscreening from 
short-range correlations, dominant at small voltages, and steric constraints of finite ion sizes, which 
prevail at large voltages. Increasing the voltage gradually suppresses overscreening in favor of the 
crowding of counterions in a condensed inner layer near the electrode. The predicted ion profiles 
and capacitance- voltage relations are consistent with recent computer simulations and experiments 
on room-temperature ionic liquids, using a correlation length of order the ion size. 

PACS numbers: 



Introduction. — The rediscovery of room tempera- 
ture ionic liquids (RTILs) as designer solvents promised 
a revolution in synthetic chemistry [l|. Thousands of 
RTILs have been synthesized with large organic cations 
and similar organic or smaller inorganic anions. Non- 
volatile and capable of withstanding up to ±4-6 V with- 
out decomposition, RTILs also hold promise as solvent- 
free electrolytes for super-capacitors, solar cells, batteries 
and electroactuators [2l4l0|. 

For such applications, it is crucial to understand the 
structure of the RTIL/electrode double layer. The clas- 
sical Gouy- Chapman- Stern (GCS) model for dilute elec- 
trolytes was used to interpret RTIL capacitance data un- 
til recently, when a mean-field theory for the crowding 
of finite-sized ions [11[ suggested bell or camel shapes 
of the differential capacitance versus voltage, decaying 
as C ^ V -1 / 2 . These were basically confirmed in sub- 
sequent experimental 12l-[l6|. theoretical [It], Eij and 
computational 2o|- 24j studies. Similar theories have also 
been developed for highly concentrated electrolytic so- 
lutions [25l427|. but none of these models accounts for 
short-range Coulomb correlations [28| . which could be 
very strong in RTIL 111 HI. As first revealed by linear 
response theories of molten salts [30[ , correlations gener- 
ally lead to over- screening [HI, where the first layer at 
the electrode delivers more counter-charge than is on the 
surface; the next layer then sees a smaller net charge of 
the opposite sign, which it again overscreens; and so-on, 
until neutrality is reached. Recent computer simulations 
of a model RTIL/electrode interface have demonstrated 
overscreening structures at low voltage, similar to exper- 
iments [29| , which are gradually overcome by the forma- 
tion of a condensed layer of counter-ions at high voltage 
[20| . as shown in Fig. [U 

In this Letter, we suggest a phenomenological theory to 
describe the interplay between over-screening and crowd- 
ing. Compared to more involved models of statistical me- 
chanics, the theory only crudely approximates discrete 
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FIG. 1: Structure of the ionic- liquid double layer (in color) 
predicted by our theory and molecular dynamics simula- 
tions [20] (Figs. 2-3 below), (a) At a moderate voltage, 

V — lOksT/e (0.26 V), the surface charge is overscreened 
by a monolayer of counterions, which is corrected by an ex- 
cess of co-ions in the second monolayer, (b) At a high voltage, 

V — lOOksT/e (2.6 V), the crowding of counter-ions extends 
across two monolayers and dominates overscreening, which 
now leads to a co-ion excess in the third monolayer. Due to 
electrostriction, the diffuse double layer (colored ions) is more 
dense than the quasi- neutral bulk liquid (white ions). 



interactions near a surface, but it is simple enough to 
be applied to dynamical problems in nanotribology, elec- 
troactuation, and porous super-capacitors. 

Theory. — We propose a Landau- Ginzburg-like func- 
tional for the total free energy (3lj : 

G = J dr{ ff + ^-|[|V<A| 2 + ^(V 2 </)) 2 ]} + j> dvqst 

(1) 

where g(c + , c_) is the enthalpy density, depending on the 
ionic concentations c±, as described below; p = e(z + c + — 
Z-cJ) is the mean charge density in the liquid volume V; 
q s is the surface charge density on a bounding metal sur- 
face S; (j) is the mean electrostatic potential, and we sub- 
tract the self energy of the electric field — ||V(/>| 2 , assum- 
ing a constant permittivity e to describe the polarizabil- 
ity of the ions. The first three terms in brackets are those 
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used in mean-field theories of ionic liquids , ionic crys- 
tals [32[ and electrolytes 25|, [33]. To go beyond that ap- 
proximation, we introduce the next allowable potential- 
gradient term, — |i^(V 2 (/>) 2 , similar to Cahn-Hilliard 
concentration-gradient expansions IS, [HI , where £ c is 
an electrostatic correlation length [31( . 

The sign of the correlation term is negative to describe 
over- screening in strongly correlated liquids: The energy 
is lowered by enhancing the curvature of </>, a measure of 
the "mean-field charge density", p = — e\/ 2 (j). For point 
charges, £ c is on the order of the Bjerrum length Ib — 
(ze) 2 /47rskT (in SI units). For RTILs with e « l(k , the 
Bjerrum length, Ib ~ 5.5 nm, is much larger than the ion 
diameter, a « 1 nm so the correlation length i c « a 
is typically at the molecular scale (29| . 

Setting 5G/5<p = for bulk and surface variations [3lj . 
we obtain a modified Poisson equation 41] and mod- 



ified electrostatic boundary condition, respectively: 
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where D is the displacement field. Due to correlations, 
the medium permittivity f, defined by D = — £V</>, is 
a linear differential operator, e = — £ 2 V 2 ), whose 
Fourier transform (valid for wavenumber \k\ <C ^^T 1 ), 
~ s(l + £ 2 k 2 ), increases with /c, as is typical for 
molten salts [36]. It is important to note that our § is 
not the complete dielectric function of the ionic liquid, 
which should diverge at small /c, as for any conducting 
medium Q . This divergence is subtracted since transla- 
tional degrees of freedom are treated explicitly via /?(</>), 
which also takes into account the nonlinear response in 
the rearrangement of ions. In our model, e approximates 
the linear dielectric response of the liquid of correlated 
ion pairs (zwitterions), which are considered to be bound 
by stronger forces, independent of the mean electric field. 

Since Poisson's equation ([2]) is now fourth-order, we 
need additional boundary conditions, similar to electro- 
dynamics with spatial dispersion [37]. Consistent with 
our bulk gradient expansion, we neglect correlations at 
the surface and apply the standard boundary condition, 
—eh-Vcj) = q s . Equation (j3]) then implies n- V(V 2 (/>) = 0, 
which requires that the mean-field charge density is "flat" 
at the surface, h ■ Vp = 0, consistent with a continuum 
model of finite-sized ions. 

Following Ref. 



the classical model 



111, we describe crowding effects via 



42]: 



k B T 



{vc+ \n(vc+) + vc- ln(vc_) 



+ [1 - v(c+ + c-)]hi[l - v(c+ + c-)]} (4) 

which is the entropy density g = —TS/v of an ideal so- 
lution of cations, anions, and holes, respectively, of min- 
imum volume v. We set v = (ir /6)a 3 /$> max = 0.83a 3 
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FIG. 2: Volt age- dependent double- layer structure predicted 
by our model, (a) Dimensionless charge density versus x/a 
for V — eV/kBT — 1, 10, 100 (solid curves), compared to the 
mean- field theory [ll| with S c = (dashed curves), (b) Di- 
mensionless cation (solid) and anion (dashed) concentrations 
and mass density (dash-dot) at high voltage, V = 100. Po- 
sition x is measured from the distance of closest approach 
and scaled to the ion diameter a = 10A. Model parameters 
7 = 0.5 (bulk/maximum density), 5 C = 10 (correlation/Debye 
length), and e — 5eo are estimated from ion profiles in simu- 
lations [U (Fig. 3 below). 



for random close packing of spheres at volume frac- 
tion &max — 0.63. More accurate expressions for g 
are available for uniform hard-sphere mixtures 25], but, 
due to the breakdown of the local-density approxima- 
tion 1281 , they over-estimate steric repulsion in the double 
layer [38]. The weaker repulsion in (J2J) actually provides 
a better first approximation for the packing entropy. 
The electrochemical potentials of the ions are then 



SG 7 
p± = -z — = k B T In 

dc± 



1 - v(c + + c_) 



± z±e(j), (5) 



and their gradients V p± produce ionic fluxes [25]. In 
equilibrium with a reference solution with = and 
volume fraction, 7 = 2vc^ = 2vc r _^^ the conditions 
H± =constant determine the Fermi-like charge density 
distribution, p{4>). In electrolytes, 7 is the volume frac- 



3 




1 2 3 4 5 6 7 



1 2 3 4 5 6 7 








lol=32.0 jiC/cm 2 







1 2 3 4 5 6 7 
i 



FIG. 3: Distributions of cations (above) and anions (below) 
sorted into monolayer bins i = 1,2,... for different surface 
charges a, as predicted by our model (solid bars) in qualitative 
agreement with simulations (Fig. 2 of Ref. |2jj, open bars). 



tion of solvated ions in the bulk [25|, [26, 33, 39]. In ionic 
liquids 7 (< 1) is the ratio of the bulk ion density to 
the maximum possible density, which characterizes their 
ability to compress [11[ . In equilibrium, we obtain a (di- 
mensionless) modified Poisson-Fermi equation, 



(1 - <5 C 2 V 2 )V 2 



sinh ( 



l + 27sinh^(^/2) 



-p($) (6) 



where x = x/\d, V = AdV, 4> = ze<p/kBT. Here, 
ze is the Debye screening length, and 
$c = ^c/^d is the dimensionless correlation length, which 
controls deviations from the mean-field theory. For 
e = lO^o and a = 10 A, the Debye length is very small, 
Ad = 1.1 A, so the ion size a becomes the relevant length 
scale [43]. If we chose 5 C = 10 to reproduce double-layer 
properties from simulations [lo| (below), then correla- 
tions are indeed at the molecular scale, £ c ~ a. 

Results. — Let us apply our model to a half space by 
solving 5l<j>""_-<j>" = p(4>) for x > subject to 0"'(O) = 0, 
and (f)(0) = V = zeV/ksT, where V is the surface poten- 
tial relative to the bulk. We solve the model analytically 
for small, moderate and large voltages [31( and compare 
with numerical solutions. 

1. Structure of the double layer. In Fig. [5] we show 
the calculated charge density (a), mass density and ion 
concentrations (b) for 7 = 1/2 and S c = 10. For a = 10 
A, T = 450K and e = 5eo, which imply i c = 0.95a, 
the model predicts molecular-scale charge-density oscil- 
lations, similar to experiments [29| and in good agree- 
ment with simulations [2o|, as shown in Fig. 03 At small 
potentials, the oscillation period and damping length are 
A ~ 2tt^/28~ c for 5 C ^> 1 (3l[, or with units restored, 
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FIG. 4: Double-layer differential capacitance Cd from our 
model (solid), simulations [20] (dashed), mean-field the- 
ory [ll[ (dash-dot), and our asymptotic scalings (inset). 



A ~ 2tt\/\d£c = 20 A= 2.0a. With increasing volt- 
age, a condensed layer of counter ions forms and expands 
into the bulk, as predicted by the mean- field theory [Tlj . 
but with the important difference that this layer over- 
screens the surface charge, leading to a second layer of ex- 
cess co- ions, which again (slightly) overscreens and trig- 
gers the same low- voltage damped charge-density oscilla- 
tions. The model also predicts non-uniform electrostric- 
tion at high voltage (Fig. Efb)) consistent with simu- 
lations (Fig. [3]): The first counterion layer attains the 
maximum density, while the next co-ion-rich layer has a 
lower density, but still larger than the bulk. 

2. Double-layer capacitance. An important property 
of the double layer is its voltage-dependent capacitance 
C(V). It has been found that excluded volume effects 
explain trends in the experimental data, but the mean- 
field theory over-estimates C, unless an empirical Stern- 
layer correction is added [20|, l21(. In Fig. [7] we show 
the double-layer capacitance versus voltage in our model, 
which is in very close agreement with simulations of Ref. 
(2o| without fitting any additional parameters. We only 
account for the extra capacitance, C s = 2s /a, in se- 
ries with the diffuse double layer, due to the distance 
of closest approach of ion centers, a/2. The value of C s 
relative to the mean- field Debye value, Cd = s/Xd, is 
C s = C S /C D = 2X D /a * 2/S c . 

At low- volt age, the model can be linearized and solved 
to find the diffuse layer capacitance, Cd [3l|, 



c d 



Cd\p V2d c + 1 
~ 6 C + 1 



for I VI < 1. 



(7) 



By extending the Composite Diffuse Layer Model of Ref. 
|26j we can also approximate C'd at moderate voltages, 
once the condensed counterion layer forms and 5 C 3> 1; 



C'd 



8 3/4 



This scaling breaks down at very large voltages when 
the condensed layer of charge grows enough to dominate 

the capacitance, yielding Cd ^ \JVl^ as m the mean- 
field theory [111, These scalings compare well with 
numerical solutions for S c ^> 1 [31] and explain why our 
model is closer to simulations than the mean-field theory 
without correlations (Fig. [7]). 

Conclusion. — In this paper we have made a first 
attempt to describe both overscreening and crowding in 
dense Coulomb liquids, such as RTILs and molten salts. 
Our simple phenomenological theory predicts that over- 
screening is pronounced at small voltages and gradually 
replaced by the formation of a condensed layer of coun- 
ter ions, followed by complete lattice saturation at very 
large voltages. Each of these three regimes is charac- 
terized by its own capacitance- voltage dependence. Our 
findings are in line with simulations and experiments, 
and they give a more complete picture of the nonlinear 
polarization of ionic liquids. 
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SUPPLEMENTARY INFORMATION 



PHENOMENOLOGICAL THEORY OF 
ELECTROSTATIC CORRELATIONS 

Let G = G e i + G c hem, where G e i is the electro- 
static energy and G c hem — J v drg is the chemical (non- 
electrostatic) part of the free energy. Suppose that G c hem 
is known, and focus on electrostatic correlation effects in 

G e l. 

The electrostatic potential, 0, is the free energy per ion 
(free charge). The electrostatic energy cost for adding a 
charge Sp in the bulk liquid volume V or Sq s on the metal 
surface S is, 



SG e l 



/ dv <j)5p - 
Jv 



\ drcj)5q s . 
Js 



(9) 



The charge is related to the displacement field D via 
Maxwell's "first" equation, 



V D 



Sp = V • 5T>. 



(10) 



The corresponding boundary condition for an ideal metal 
surface (where D = 0) is, 



[n • D] = n • D = —q s 



Sq s 



D. 



(11) 



Substituting these expressions into ([9]) and using Gauss' 
theorem, along with the definition of the electric field, 
E = — V(/>, we recover the standard electrostatic free en- 
ergy equation [? ] , 



SG e 



drE-£D. 



(12) 



In the linear response regime (for small external elec- 
tric fields), we have 



D = £E, 



where e is a linear operator, whose Fourier transform s(k) 
encodes how the permittivity depends on the wavelength 
2tt jk of the /c-Fourier component of the field, due to dis- 
crete ion-ion correlations, as well as any non-local dielec- 
tric response of the ions, such as exponentially decaying 
Debye correlations in ionic plasma, as well as correlations 
in polarization flucutations due to any other molecules if 
they are present in the liquid. We can then integrate ([T2]) 
over ST) through a charging process that creates all the 
charges in the bulk and surface from zero to obtain 



G e . 



drE D. 



(14) 



For a given distribution of charges p and q s , with as- 
sociated displacement field D, the physical electric field 



E is the one that minimizes G e i-> subject to the con- 
straint of satisfying Maxwell's equations ([T0]) - ([TT]) . Since 
E = -V0 to enforce V x E = 0, we can minimize G e i 
with respect to variations in (/>, using Lagrange multipli- 
ers Ai and A2 to enforce the constraints, 



G el [<t>] 



dr 



E • D + Ai (p - V • D) 



+ j> dr s A2 (q s + n ■ D) 



(15) 



To calculate the extremum, we use the Frechet functional 
derivative: 

SGel v Gel [<t> + e</> £ e ] - Gel [4>] ~n 

= lim 7 (I 6 ) 

S(j) e^O 600 

where S(j) € = (j) 5 e (r, r') is a localized perturbation of the 
potential (with compact support), which tends either to 
a 3D delta function in the liquid (r G V) or to a 2D delta 
function on the surface (r G S) as e — >> 0, and (j)o is an 
arbitrary potential scale for dimensional consistency. By 
setting 5G e i /5(f) = for both surface and bulk variations, 
we find Ai = A2 = </>. Finally, using vector identities, we 
arrive at a general functional for the electrostatic energy, 



G el [<j>] 



dr 



dr s q s (j) (17) 



to be minimized with respect to </>, once we know the 
relationship between D and E = — V</>. 

To model the field energy in an ionic liquid, we assume 
linear dielectric response of the molecules with constant 
permittivity e plus a non-local contribution for ion-ion 
correlations. Here, the permittivity e describes the elec- 
tronic polarizability of the ions. 



9 field --V0-D 



(13) where 



E(r) 2 + f dv'K(v,v')p(v)p(v') 
Jv 

(18) 



eV • E 



(19) 



is the "mean-field charge density", which would pro- 
duce the electric field in the dielectric medium without 
accounting for ion-ion correlations. Suppose that the 
non-local kernel K(r, r') decays over a length scale £ c , 
bounded below by the finite ion size a and above by the 
Bjerrum length £5, which sets the scale for electrostatic 
correlations among point charges. For charge variations 
over scales larger than £ c (corresponding to small pertur- 
bation wavenumbers, £ c \k\ <C 1), we obtain a gradient 
expansion for the non-local term 



9 field 



m 2 + E 



71=0 



tT 



-V n p 



(20) 
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where a n are dimensionless coefficients, which implies 



G e 



dr 



|V0| 2 + 5>»- 2 (*T 1 v n < 



dr.. 



(21) 



to position x from the surface versus x. The integrated 
charge is then normalized by the total charge in the dou- 
ble layer, as in Ref. [2o|. This graph provides a quanti- 
tative characterization of the strength of over-screening 
in the first layer. 



Equation (1) in the main text results from the first term 
in the gradient expansion of the non-local electrostatic 
energy with the choice ao = 1 (after suitably rescaling 
£ c ), where the overall negative sign of this term is chosen 
to promote over-screening. 

By settting SG e i / 5(j) = for bulk and surface perturba- 
tions in (|2T|) . we recover Maxwell's equations ([T0]) -([TT ]) . 
with D = £E, where the permittivity operator has the 
following gradient expansion, 



i = e l-^ a „-^V 2n 



(22) 



and corresponding small- k expansion of the Fourier 
transform, 



i(k) 



l + ^an-iC-l)"" 1 ^) 2 " (23) 

71=1 

[1 + a (£ c k) 2 ] (24) 



which grows with k at small wavenumbers in the case 
where correlations promote over screening, ao > 0. Note 
that it is well known that such an expansion only holds 
at small k. At larger fc, e(k) diverges, becomes negative 
on the other side of the singularity, then diverges again 
to — oo at another point, and becomes positive after the 
second divergence; see Refs. (iol. l36|. 



CHARGE PROFILES AND VERIFICATION OF 
OVER-SCREENING 



LOW- VOLTAGE ANALYTICAL SOLUTION 

The solution to the equations must be calculated nu- 
merically. However, at low voltage we can obtain an an- 
alytical solution. At low voltage, we have the approxi- 
mation that, 



' dx 4 



(25) 



The analytical solution to this equation depends on 
whether S c is greater than, equal to, or less that \. Since 
S c is presumed large in case of ionic liquids, we present 
the analytical solution for S c > |, 



(j){x) = Ve~ klX (cos(k 2 x) + A sm(k 2 x)) , (26) 



where 
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The total charge in the diffuse double layer can be eval- 
uated from 
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The diffuse layer capacitance in the limit when S c is large 
is approximately Cd = \/2/S c . The diffuse layer capaci- 
tance is less than the classical theory without correlations 
and decreases with the square root of S c . 



In the main text we show charge density profiles for 
a specific set of parameters. The charge density profiles 
in the text are presented in spatial coordinates scaled by 
the ion size. However, the natural length scale for the 
dimensionless problem is the Debye length. The solu- 
tions to the equation in dimensionless form depend upon 
the applied voltage, the correlation length scale 5 C1 and 
the volume fraction 7. In Fig. [5] we show the charge 
density as a function of distance (normalized by the De- 
bye length) for increasing values of S c . To convert these 
ion profiles to dimensional form, the x-axis need only be 
scaled by the value of /a as given by the physical pa- 
rameters of the problem. Fig. [5] shows that the strength 
of the over-screening is a strong function of S c . 

To prove that our simple continuum model predicts 
over-screening by the first, condensed layer of counteri- 
ons, in Fig. [6] we plot the integrated charge density up 



HIGH- VOLTAGE COMPOSITE 
APPROXIMATION 

In ionic liquids, the parameter 7 is on the order of unity 
and excluded volume effects are significant. At voltages 
beyond the linear response we find that a condensed layer 
of counter-ions forms near the wall. In this condensed 
region close the wall we could solve, 

. 2 d 4 ^ d 2 <A _ _ 

)c dx^~dx^J _/W 

where we assume that the charge density is a constant, 
Pmax, and has reached the maximum value defined by 
the value of 7; i.e. p m ax = I/7 if we apply a negative 
voltage. In order to further simplify the approximation, 
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FIG. 5: Charge density profiles for 7 = 1 and 5 C = 2, 10, and 50, as indicated. The solid curves are calculated from our model 
while the dashed curves are calculated for the case where S c — and there are no correlations. Solutions are shown for applied 
voltages of V — 1, 10 and 100 measured in units of ksT/e ~ 25 mV. Similar results are obtained for different values of the 
volume fraction, 7. 




FIG. 6: Cumulative charge density profiles as a function of distance from the electrode. The charge is integrated cumulatively 
in space and normalized by the total double layer charge; namely q(x) = J* p(x)dx/ J °° p(x)dx. Solutions are shown for 7 = 1 
and 5 C = 2, 10 and 50, as indicated. Applied voltages are V = 1, 10 and 100 measured in units of ksT/e « 25 mV. These 
cumulative profiles clearly show over-screening. 



we can assume that in ionic liquids, S c is typically large 
and we solve as an approximation, 

in the wall region. 

The general solution for the potential in the wall region 
becomes a fourth order polynomial. Using the boundary 
conditions that we have fixed voltage V at x = 0, along 
with d 3 (/)/dx 3 = at the wall, our solution for the po- 
tential has a simple form, 

^) = fef x A + Bx 2 + Cx + V 
24£;? 

This polynomial solution which is valid near the wall 
can be matched to the low voltage solution provided in 
the previous section. Ensuring continuity of the poten- 
tial, the charge density, and all the derivatives allows us 
to solve for the unknown constants of integration. The 
resulting analysis yields a quartic equation for the size of 
the condensed layer. Once the size of the condensed layer 
is known, all the constants for the matching are easy to 



obtain. The approximate composite model was found to 
match the full numerical simulation as long as the voltage 
was low enough that a second condensed layer of opposite 
charge did not begin to form (see Fig. lc at V = 100). 

While this analysis may be useful, the resulting quartic 
equation does not provide a simple form for the double 
layer capacitance. A much simpler form of this compos- 
ite solution emerges if we make the additional approx- 
imation that all of the voltage drop occurs across the 
condensed layer and, after the condensed layer, the po- 
tential and all its derivatives go to zero. While this is not 
true, it is found from the numerical solutions to provide 
a reasonable prediction of the capacitance at "interme- 
diate voltages", where a condensed layer forms, but not 
so high that a second condensed layer of opposite charge 
forms due to over-screening. While this assumption does 
not yield complete charge density profiles, it does provide 
a useful approximation for the potential in the condensed 
layer and thus the capacitance. 

Returning to the general solution and using the 
simpified boundary conditions, (j>(x = L) = and 
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FIG. 7: (a) Capacitance of the diffuse double layer, nor- 
malized by the Debye value, as a function of voltage for 
5 C = 0, 2, 10, and 50 while holding 7 = 1. Numerical 
solutions (solid curve) are compared with the low-voltage 
and moderate voltage approximations in Eqs. (7) and (8) 
of the main text (dashed lines) and the high-voltage scaling, 

Cd ~ yJl/^V (dash-dot line). We show only the diffuse layer 
capacitance to highlight the different scaling laws, (b) Total 
capacitance of the double layer for 5 C = 10 and 7 = 1; the 
inner layer is included in series with diffuse layer. 
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= 0, we find, 
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Setting the second derivative to zero at x 
size of the condensed layer, 



L= (Vj5 2 c 8)*. 



L yields the 



Solving for the total charge q = -S 

x x—0 



we obtain a 

simple approximation for the diffuse layer differential ca- 
pacitance (Cd = dQ/dV), 
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(28) 



The scaling presented above is essentially valid at large 
S c and moderately large voltages. At high voltage, cor- 
relations become irrelevant, since "crowding beats over- 
screening", and the capacitance is determined by the 
excluded volume effects only and has a scaling Cd ~ 

as previously discovered. 

The transitions between these three regimes are evi- 
dent in Fig. [3 where we compare the numerical solution 
for the capacitance to the simple scaling laws derived 
above. The range of validity of the intermediate volt- 
age expression simply comes from the intersection of the 
three regimes. In Fig. [7^ we show the diffuse layer ca- 
pacitance only to clearly show the results of the simple 
scaling laws. At 5 C = 10 we find a short transition regime 
where the V~ X ^ A scaling appears, and for S c > 50 the 
scaling is valid over a wide range. All the capacitance 
curves appear to converge to the S c = solution at high 
voltages, as expected from our analysis. When S c is large, 
however, the voltages where the capacitance curves con- 
verge are quite extreme, so this limiting behavior may 
have limited applicability. Perhaps it could be used to 
validate simulations. 

In Fig. [7)3 we show the total capacitance (inner layer 
and diffuse layer in series) for S c — 10 and 7 = 1. Here 
we find that the simple scaling laws are useful for under- 
standing the capacitance computed from the numerical 
solution of our model. 




1 2,3 4 5 

x/a 



